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ABSTRACT 

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive 
method for measuring the build-up of cosmological structure, which depends both on the ex- 
pansion rate of the Universe and our theory of gravity. The statistical precision with which 
redshift space distortions can now be measured demands better control of our theoretical sys- 
tematic errors. While many recent studies focus on understanding dark matter clustering in 
redshift space, galaxies occupy special places in the universe: dark matter halos. In our de- 
tailed study of halo clustering and velocity statistics in 67.5/z~ 3 Gpc 3 of A^-body simulations, 
we uncover a complex dependence of redshift space clustering on halo bias. We identify two 
distinct corrections which affect the halo redshift space correlation function on quasilinear 
scales (~ 30 - 80 /r'Mpc): the non-linear mapping between real and redshift space positions, 
and the non-linear suppression of power in the velocity divergence field. We model the first 
non-perturbatively using the scale-dependent Gaussian streaming model, which we show is 
accurate at the < 0.5 (2) per cent level in transforming real space clustering and velocity statis- 
tics into redshift space on scales s > 10 (s > 25)/z~ 1 Mpc for the monopole (quadrupole) halo 
correlation functions. The dominant correction to the Kaiser limit in this model scales like b 3 . 
We use standard perturbation theory to predict the real space pairwise halo velocity statistics. 
Our fully analytic model is accurate at the 2 per cent level only on scales s > 40/i~'Mpc for 
the range of halo masses we studied (with b — 1.4 — 2.8). We find that recent models of halo 
redshift space clustering that neglect the corrections from the bispectrum and higher order 
terms from the non-linear real-to-redshift space mapping will not have the accuracy required 
for current and future observational analyses. Finally, we note that our simulation results con- 
firm the essential but non-trivial assumption that on large scales, the bias inferred from real 
space clustering of halos is the same one that determines their pairwise infall velocity ampli- 
tude at the per cent level. 

Key words: cosmology: large-scale structure of Universe, cosmological parameters, galaxies: 
haloes, statistics 



1 INTRODUCTION 

The growth of large-scale structure, as revealed in the clustering of 
galaxies observed in large redshift surveys, has historically been 
one of our most important cosmological probes. This growth is 
driven by a competition between gravitational attraction and the 
expansion of space-time, allowing us to test our model of grav- 
ity and the expansion history of the Universe. Despite the fact that 
galaxy light doesn't faithfully trace the mass, even on large scales, 
galaxies are expected to act nearly as test particles within the cos- 
mological matter flow. Thus the motions of galaxies carry an im- 
print of the rate of growth of large-scale structure and allows us to 
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both probe dark energy and test General Relativity (see e.g. |Jain & 
| Zhang 2008 ; Nes seris & Perivolaropoulos|2 008 , Song & Koyama 
2009[|Song & Percival|2009||Percival & •WMte|2008||McDonald"& 
|seljak|2008||White, Song & Percival|2009| for recent studies). 

This measurement of the growth of structure relies on redshift- 
space distortions seen in galaxy surveys (Kaiser 1987). Even 
though we expect the clustering of galaxies in real space to have 
no preferred direction, galaxy maps produced by estimating dis- 
tances from redshifts obtained in spectroscopic surveys reveal 
an anisotropic galaxy distribution. This anisotropy arises because 
galaxy recession velocities, from which distances are inferred, in- 
clude components from both the Hubble flow and peculiar veloc- 
ities driven by the clustering of matter (see Hamilton 1998, for a 
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review). Measurements of the anisotropies allow constraints to be 
placed on the rate of growth of clustering. 

On large scales, where linear perturbation theory is valid, 
it is natural to work in a Fourier basis because the symmetries 
of the background solution imply that k-modes evolve indepen- 
dently. On smaller scales, and especially once survey non-idealities 
and fingers-of-god become important, the choice is not so clear. 
Because the velocity field departs from its linear theory predic- 
tion on extremely large scales (k < 0.03/iMpc~'), models be- 
yond linear theory must be used to extract cosmological informa- 
tion from redshift surveys. This has been long recognised and a 
variety of methods have been attempted to model the distortions. 
Several recent studies of redshift space distortions have provided 
non-linear descriptions of the matter density field in Fourier space 
l|Taruya, Nishimichi & Saito||2010[ |Iennings, Baug h7"& Pascoh| 
|2011a| > which agree well with direct N-body calculations of the 
effect. However, we do not generally observe the matter density 
but rather tracers which tend to live in dark matter halos. This in- 
troduces further effects which must be carefully modelled if we 
are to achieve the desired accuracy. In this paper we find a strong 
dependence on halo bias in the shape of the redshift space correla- 
tion function, indicating the need for more sophisticated theoretical 
models (see also |Taruya, S aito & Nishimichi 2010, for a Fourier- 
space approach). We trace this strong bias dependence primarily 
to the non-linear mapping between real and redshift space. While 
more slowly varying with bias, non-linear evolution of the pairwise 
velocity distribution also substantially changes the redshift space 
clustering in comparison to linear theory. 



[Tinker, Weinberg, & Zheng (2006) and |Tinker|p007) , build- 
ing on the work of Hatton & Colej|l999), combined the streaming 
and halo models to describe the redshift space correlation function 
on scales of r < 20 fr~'Mpc. While we take a similar approach here, 
we focus on larger scales and only on halos and ignore the contri- 
butions from satellite galaxies for now. This greatly simplifies our 
modelling compared with Tinker (20071; however, as we will see, 
there is still a rich phenomenology compared to dark matter clus- 
tering. 

The outline of the paper is as follows. We first review linear 
and quasilinear descriptions of redshift space distortions in both 
Fourier and configuration space, and introduce the scale-dependent 
Gaussian streaming model that we study in detail in this paper. In 
Section [5] we present a simplified Fisher matrix calculation of the 
expected constraint on the peculiar velocity field from the ongoing 
Baryon Oscillation Spectroscopic Survey (BOSS; Schlegel, White, 
|& Eise nstein 2009, Eisenst ein, et al.|2011) . This calculation sets 
the target for the accuracy of our model. In Section|4]we compare 
the streaming model ansatz as a transformation between real space 
clustering and velocity statistics and clustering in redshift space 
for halo mass bins, finding very good agreement. Then Section [5] 
examines the ability of perturbation theory to describe the four in- 
gredients of the streaming model ansatz for tracers of the matter 
density field: the real space correlation function, the mean tracer 
pairwise velocity, and the tracer velocity dispersions along and per- 
pendicular to the line-of-sight (LOS). We show in Section [6] that 
using our perturbation theory description as input into the scale- 
dependent Gaussian streaming model, we have an analytic model 
accurate at the 2 per cent level on scales s > 40 hT 1 Mpc. Sections]?] 
and [8] identify the dominant non-linear terms, elucidating how our 
model depends on redshift and halo bias. 



2 MODELLING REDSHIFT SPACE DISTORTIONS 

We begin by reviewing the effect of redshift space distortions in 
linear theory both in Fourier space and configuration space. While 
Fourier space is usually preferred for theoretical investigations, 
configuration space is simpler when dealing with wide-angle ef- 
fects (e.g. |Papai & Szapudi 2008 1 and often preferred in observa- 
tional analyses. 



2.1 Linear theory: Fourier space 

The redshift-space position s of a galaxy differs from its real-space 
position r due to its peculiar velocity, 



s = x + v z (x) z, 



(1) 



where v-(x) = u z (x)/(aH) is the line-of-sight (LOS) component of 
the galaxy velocity (assumed non-relativistic) in units of the Hubble 
velocity, and we have taken the LOS to be the z-axis. We shall adopt 
the "plane-parallel" approximation, so this direction is fixed for all 
tracers (halos, galaxies, etc.). 

The galaxy over-density field in redshift-space can be obtained 
by imposing mass conservation, (1 + 6 s g )d 3 s = (1 + d g )d 3 r. For 
a uniform, z-independent mean galaxy density, the exact Jacobian 
for the real-space to redshift-space transformation is 

d 3 s / v.\ 2 / dv.\ 

!l + 7 (2) 



cPr 



In the limit where we are looking at scales much smaller than the 
mean distance to the pair, vjz is small and it is only the second 
term that is important (Kaiser 1987] but see |Papai & Szapudi|2008| 
|Shaw & Lewis|2008| >, 

-1 



1 + 



dv z 
~dz 



(3) 



If we assume an irrotational velocity field we can write v ; = 
-3/3zV~ 2 6>, where = -V • v, and V~ 2 is the inverse Laplacian 
operator. In Fourier space, (d/3z) 2 V~ 2 = (k-Jk) 2 = fi 2 , where is 
the cosine of the LOS angle, so we have that 



(4) 



to linear order. Often it is further assumed that the velocity field 
comes from linear perturbation theory. Then 



0(k) = /cWss(fc), 



(5) 



where / = dlnD/dlna * Q~ 6 {Peebles 1980 1 



For a population of galaxies, which we denote with a subscript 
g, the linear, redshift-space power spectrum is then proportional to 
the linear, real-space matter power spectrum, P' m (k) 



P s g (k,fi) = (b + fff P r Jk) = b 2 (l +Pti 2 ) 2 P r m (k). 



(6) 



2.2 Linear theory: configuration space 



Fisher (1995 ) first showed in detail the relation between the Kaiser 
formula in Fourier space, Eq. [6] and the redshift space correlation 
function; we rely heavily on that work for this section. In linear 
theory, the correlation between 6(x) and v(x') gives rise to a mean 
infall vi2(r) between pairs of matter tracers. The velocity disper- 
sion along the LOS, cr 2 2 (r) = ((v z (x) - v z (x') 2 ^, depends both on 
scale and the orientation of the pair separation vector with respect 
to the LOS. These scale dependencies give rise to linear redshift- 
space distortions in configuration space, apparent in Fig.[TJas the 
squashing of contours along the Z axis. 
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Figure 1. Contours of constant f as a function of separation perpendicular 
(R) and parallel (Z) to the line-of-sight (LOS) from the mock galaxy cat- 
alogues in White et al.|(20TT) . The magnitude of the squashing along the 
Z axis depends on the amplitude of the peculiar velocity field, which is re- 
lated to the over-density field on large scales (Eq.[5). The stretching of the 
contours near R=() results from small-scale velocity dispersions and is often 
referred to as the 'fingers-of-god.' 



We reproduce the relevant linear theory velocity predictions 
here, but generalised to linearly biased tracers with bias b: 



Vi2« = v 12 (r)r = -r^ J dk k P' m (k)j ,(Jfcr) 
(v,(r' + r)Vj(r')} = T ± (r)^ + [>P„(r) - ^A^Wj 

*x(r)=^ f dkp;„(k) J - l(kn 



jo(kr) - 



kr 

Vim 



kr 



a\ z {r tl i 2 ) = 2 [o? - M 2x f\\(r) - (1 - SWAr)] 



(7) 
(8) 
(9) 

(10) 
(11) 



where \i = 1 ■ r denotes the angle between the LOS and the pair 
separation vector and a 2 is the one-dimensional velocity disper- 
sion, j (v(x) • v(x)). Our Fourier convention is 



J (S^*'V"(k); /(k) = J d 3 xe-**f(k). (12) 



|Fisher| ( fl995| > derived the linear theory redshift space distor- 
tion limit in configuration space under the assumption that the den- 
sity and velocity fields are Gaussian, which amounts to evaluating 
the following expression: 

1 + e„(r<r, r ^ = \f dy a+6i)(l+S 2 ) S D (r„ - y + v, - v 2 ) 

'(13) 

where Sip and v 1/2 are the density and LOS velocity at two points 
with real space separation y along the LOS and r a perpendicular to 
the LOS, and 6 D is the Dirac delta function ensuring that the pair is 
mapped to redshift space separation r n . We can then re-express 6° 



\+?Jr*,r K ) = [ dy (1 +tf 1 )(l + fc) 



dK 
2n 



(14) 



and compute the expectation value, assuming Gaussian statistics 



(eq. 20 of |Fisher|1995) : 

l+?f(rr,r w ) = J 



dy 



exp 



(r» - yf 



l+f t (r) + 



y (r x -y)vn(r) 1 y 2 v l2 (r) 
4 r 2 <r 2 2 (y) 



1 - 



(15) 



If we expand Eq.[T5]to linear order, the redshift space correlation 
function is equivalent to 



e t (r„r x ) =f f (s)- j\>nir) 



+ - — [o- 2 l2 (y) 



2 dy 2 



(16) 



where ^ r g (s) = b 2 i; r m {s) is the real space linear galaxy correlation 
function evaluated at redshift space separation s 2 = r 2 . + r\ and 
r 2 = r 2 a + y 2 is the real space separation of the pair. We expand 
the last two terms to elucidate the contribution from different halo 
velocity statistics as a function of angle and separation. The terms 
proportional to bf in Eq. [6] arise from the pairwise infall and its 
derivative as a function of separation r (denoted with ' through- 
out), and terms proportional to f 2 arise from the pairwise velocity 
dispersion and its first and second derivatives: 



I l Vn(r) r 



Vn(.r) 



(17) 



2dy 2 



[o- 2 2 0>,/v)]| _ 



(2- 



10// + 8//) , " ; 



2 4 V 2 

+ 5fi 4 )-^— + (-1 + 6fi 2 



, Mr) 

5yU ) 

r r 



(18) 



Since the theoretical prediction for the redshift space correlation 
function depends on first and second derivatives of velocity statis- 
tics, scale-dependent errors on a theoretical model for these func- 
tions can translate into large errors on £ s . Eq.[T8]also demonstrates 
that at linear order, any constant, isotropic velocity dispersion does 
not alter the redshift space correlation function, since only the dif- 
ference x ¥±(r) - YyM enters, along with derivatives of those func- 
tions. 



2.3 Moments 

It is standard practise to expand the dependence of P(k) or <f(s) on 
the LOS angle in Legendre polynomials, which we here write as L{ 
to avoid confusion with the power spectrum moments: 



P(k,fi k ) = Y j P e (k)L e (p k ). 



The power spectrum moments are given by 
2£+l 



Pt(k) = 



I 



dfi P(k,n k )L e (Mk)- 



(19) 



(20) 



In linear theory, only £ = 0, 2, and 4 are non-zero. The power spec- 
trum is particularly simple because the ^-dependencies of the Leg- 
endre moments P[(k) are identical, and the relative amplitudes de- 
pend only on the sample bias b and the rate of structure growth 
/■ 



(Po(k)) 


1 




= PUk) 


W)J 


V 



(21) 
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2.4 Beyond linear theory 
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Figure 2. Comparison of A'-body simulation redshift space clustering of 
halo samples in Table[T]with the biased Lagrangian perturbation theory of 
Matsubara (2008b). The high (circles) and low (triangles) halo mass bins 
are shifted by 10 per cent for clarity compared with the HOD subsample 
(diamonds) in the upper left panel. The predictions of Matsubara 1 2008b) 
are shown as solid black curves; higher bias values have larger deviations 
from linear theory expectations (dashed curves) in LPT. While M atsubara] 
1 2008b 1 provides an excellent fit to the simulation results for fo and £j on 
BAO scales (bottom right panel), the quasilinear scale description is not 
accurate. 



Fourier transforming to configuration space gives (e.g. |Cole et aT| 
[r^|HamTiton|1998| > 



(k)j c (ks) (22) 



where j e is the spherical Bessel function of order £ and A 2 , = 
k 3 P ( (k)/(2n 2 ) is the moment of the dimensionless power spectrum 
or power per Ink. While the relative amplitudes of the Legendre 
moments are also given by Eq.[2T| in configuration space the 
moments depend differently on s (see solid lines in Fig.[4](. For ex- 
ample, using the recurrence relations among the spherical Bessel 
functions we can write f 2 as the average value of £ up to s minus 
the value at s. 

Non-linear effects, and especially fingers-of-god on small 
scales, can cause the £ > 4 moments to be non-zero. The fingers-of- 
god are evident in Fig.[T|as the stretching of contours localised near 
R = 0. The presence of strong fingers-of-god also drives £ 2 posi- 
tive, whereas super-cluster infall on larger scales generally predicts 
fa < 0. Generally the monopole (C = 0) and quadrupole ({ = 2) 
dominate in terms of signal-to-noise, and the ratio P2/P0 ( or an 
analogous combination of £0 and £2) encodes information about 
P — f lb (but see Berlind, Narayanan & Weinberg 2001 , for a dis- 
cussion of issues with this approach). 



As is so often the case, our strongest statistical constraints come 
from the smaller scales where we have a large number of indepen- 
dent samples within our survey. These are also the scales for which 
the simple linear theory described in the last section is the least ap- 
plicable. A standard approach to include small-scale non-linearities 
manifest as fingers-of-god is to use a streaming model, where lin- 
ear theory is spliced together with an approximation for random 
motion of particles in collapsed objects (e.g. Peacock 19921. This 
arises from ignoring the scale-dependence of the mapping between 
real and redshift space separations and assuming an isotropic veloc- 
ity dispersion and amounts to convolving the linear theory result 
with a LOS smearing or multiplying the power spectrum by the 
Fourier transform of the small-scale velocity PDF. Peacock (1992 ) 
modelled the small-scale velocity field as an incoherent Gaussian 
scatter which amounts to 



(b + fffrjfie 



(23) 



A more realistic distribution for the pairwise velocities is exponen- 
tial (Peebles 1976; Davis & Peebles 1983), because galaxies popu- 
late halos of a wide range of masses and velocity dispersions ( Sheth 
[T996l|Diaferio & Geller 1 1996] | White |200 1 [| Selj ak|200 1 ^ . 

While this extension from linear theory can improve agree- 
ment with observations and A'-body simulation results, this sim- 
ple approach ignores the well-known fact that the velocity diver- 
gence field 9 deviates from linear th eory on extremely large scales 



(its 0.03 fcMpc -1 ); see e.g., fig. 8 of|Carlson, White, & Padmanab 



Scoccimarro| ( |2004 1 therefore 



|han| ( [2009| > for a recent comparison, 
proposed a simple extension of linear theory for the matter power 
spectrum in redshift space: 



Di.PT 1 1 \ , 1< 2nr,PT , A 4 or.PT 

\P 6S (k) + 2fn P' +f tiP m 



-k 2 ti 2 crl 



(24) 



where the power spectra on the right hand side are evaluated in 
perturbation theory in real space, and a\ is the linear velocity dis- 
persion. Jennings, Baugh, & Pascoli ( 20 1 1 a bi) have extended this 
model for the redshift space matter power spectrum to different 
dark energy scenarios by fitting a relation between the three non- 
linear power spectra in Eq.|24]to A'-body simulation outputs, but 
treat a\ as a free parameter (in a similar fashion to earlier work by 
IHatton & Cole|l999> . 

However, the relation between the real and redshift space 
statistics depends on non-trivial correlations at high orders, which 
Eq. [24] ignores; many terms from the standard perturbation the- 
ory expression for the redshift space power spectrum have been 
dropped (see Section IIIB and Appendix B of Matsubara 2008a). 
However, it is well-known that the standard perturbation theory ex- 
pression for the redshift space clustering of matter is inaccurate 
(e.g., Scoccimarro, Couchman, & Frieman ( 1999); see also Mat- 
subara (2008a) and Blake, et al. (20 ll I for recent comparisons to 
A'-body simulations and observed galaxy power spectra, respec- 
tively). The importance of higher-order terms can be seen most eas- 
ily from the expression for the correlation function, Eq.[l4] where 
we see that going beyond linear theory involves many moments of 
5 and v.. Taruya, Nishimichi & Saito (2010) showed that the bis- 
pectrum terms have oscillatory features which affect the BAO fea- 
ture, and used them to obtain a better fit to A'-body simulation data 
for the matter redshift space power spectrum. However, |Taruya,| 
Nishimichi & Saito (2010) still needed to introduce a Gaussian 
damping term with free parameter crj in order to fit the A'-body 
data well. |Tang, Kayo, & Ta kada (20ll ) take a different approach, 
and try to reconstruct P r 55 (k) and P r sg (k) from P s (k) for both matter 
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sample log(M) range 



bit,, 



b LPT n (h 'Mpc) 



high 
low 
HOD 



>13.387 
12.484 - 12.784 



2.67 
1.41 
1.84 



2.79 
1.43 
1.90 



7.55 x 1(T 5 
4.04 x 10~ 4 
3.25 x 10~ 4 



Table 1. Halo samples from A'-body simulations. The first two samples 
are defined by sharp cuts in FoF halo mass (linking length 0.168), while 
the third sample contains all halos that host at least one galaxy in the mock 
catalogues of White et al. 1 20lTJ. This is equivalent to a catalogue of all 
mock central galaxies and no satellite galaxies. The bias values S;,„ and 
1>LPT we list are derived by fitting the real space halo correlation function 
to the linear theory and Lagrangian perturbation theory predictions for £(r), 
including scales r > 30A~'Mpc. b\ in assumes a scale-independent linear 
bias, while bipj is fit including the second order bias implied from the 
peak background split (see Matsubara 1 2008b I for details). 




-20 -15 



-10 -5 

Az Mpc/h 



and halos in A'-body simulations, assuming a relation like Eq. [24] 
but with a more general multiplicative damping term. Their proce- 
dure does not work for halos unless they account for a higher order 
term in perturbation theory which scales like b 2 B sst) . 

Going beyond linear theory perturbatively introduces extra 
"mode coupling" terms, some of which can be resummed into a 
^/-dependent damping (Bharadwa j|1996||Eisenstein, Seo <Sf White 
2007| |Crocce & Scoccimarro||2008[ |Matsubara||2008a|b[ |Taruya,| 
Nishimichi & Saito 2010[>. For example, in the formalism of Mat- 
subara ( 2008a b) the damping is exponential with an angular de- 



pendence (k + fkfii.) 2 or [1 + /(/ + 2)yU 2 ] k 2 and the mode-coupling 
terms go up to yu 8 . However the range of validity of perturbation 
theory is limited, and it clearly cannot be extended to the dynam- 
ics of collapsed objects, so it is not clear how much guidance these 
forms provide. 

Another issue of utmost importance in the application to 
galaxy surveys, but that was not addressed in most of these re- 
cent perturbation theory studies, is the inclusion of bias. |Heav-| 
ens, Matarrese, & Verde ( 1998) computed the redshift space power 
spectrum for biased tracers in standard perturbation theory. This ap- 
proach is complicated by the need to introduce a smoothing scale 
by which to define local Eulerian biasing, and |Roth & Porciani| 
< |20 1 1 ) > found that this biasing scheme is not very accurate com- 
pared with halos in A'-body simulations. Finally, standard pertur- 
bation theory cannot be used to predict the behaviour of the corre- 
lation function because P^ p '° of, (fc) diverges for large k, which con- 
tribute to £(r) even for large r (see discussion in Section IIIB of 
|Matsubara|2"008al >. 

Recently Matsubara ( 2008a b I have addressed both of these 
shortcomings of SPT by developing a Lagrangian perturbation the- 
ory (LPT) description including non-linear local Lagrangian bias in 
redshift space which does not require the introduction of a smooth- 
ing scale, and is well-behaved at large k for £(r) to be computed. 
In Fig. [2] we compare the prediction from this theory to our A'- 
body simulation results for halo clustering for three different halo 
masses given in Table[T] The bottom right panel shows good agree- 
ment for both the monopole and quadrupole on BAO scales (see 
also |Padmanabhan, White & Cohn|2009| |Noh, White & Padman-] 
|abhan|2 009 l. The upper left panel shows reasonable agreement for 
£o between the theory and simulations. However, the theory pre- 
dictions for £2 and £4 do not fit the quasilinear scales of interest in 
this paper (~ 30 - 80 h Mpc), and so a more accurate theory must 
be developed in order to infer the peculiar velocity field amplitude 
from halo clustering in redshift space. 



Figure 3. The radial pairwise halo velocity probability distribution func- 
tion for the HOD halo subsample in Tablelllfor pair separations 30 tr l Mpc 
< r < 31.5 /i~'Mpc (solid curves) compared with a Gaussian distribu- 
tion (dashed curves) with the same mean (-1.6/i~'Mpc) and variance 
18 (/r'Mpc) 2 . The normalisation is arbitrary. While the halo PDF has clear 
skewness and kurtosis, the PDF for dark matter particles in our simulation 
has 30 per cent larger variance and exponential tails. The mean infall be- 
tween dark matter particles (-0.9A~'Mpc) is smaller than the more highly 
biased halo samples. 



2.5 A non-perturbative real-to-redshift space mapping: the 
scale-dependent Gaussian streaming model 

As perturbation theory does a good job of describing P ee and P m 
on the relevant scales ( Carlson, White, & Padmanabhan 2009 , see 
also our section |5.1fr we conjecture that the failure of perturba- 
tion theory descriptions of the redshift space power spectra can 
be traced to the inaccuracy of a perturbative description of the 
real-to-redshift space mapping ([Scoccimarro, Couch man, & Frie^| 
man 1999). Fisher ( 1995) derived the exact result for the redshift 
space correlation function in the case where both density and ve- 
locity fields are Gaussian and related to one another as in linear 
theory (our Eq.|15}. Thereafter, several authors (Bharadwaj 2001; 
Scoccimarro 2004 , Shaw & Lewis 2008 1 showed that there are sig- 
nificant differences between the exact result and the Kaiser limit, 
Eq. [6] which can be traced to additional assumptions inherent in 
the Kaiser derivation. However, we found that this exact mapping 
for Gaussian fields does not improve agreement with A'-body sim- 
ulation results for halo clustering in redshift space. 

A principle object of interest in the study of redshift space 
clustering statistics is the pairwise velocity probability distribution 
function (PDF) P(v z ,r), i.e., the probability that a pair with real 
space separation r has relative LOS velocity v.. Even for the ex- 
act result in the Gaussian case (Eq. |1S) , the corresponding pair- 
wise velocity PDF is non-Gaussian; however, near its peak it can 
be approximated by a Gaussian centred on yuv 12 (r) ( |Scoccimarro| 
2004). By resumming the mean infall velocity term (oc jxvu{r)) 
from Eq. [T5] into the exponential, we recover a streaming model 
expression that agrees to linear order with the Gaussian case, and 
for which the corresponding pairwise velocity PDF's mean and dis- 
persion are correct: 



l+^OV,^): 



Vr„-y- l iv n (r)Yl2o-Ur,ii)_ 



dy 



(25) 
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Figure 4. The solid curves show the linear theory Kaiser formula predictions for the Legendre moments of the redshift space correlation function, fo,2.4(s) 
given in Eq. |22| where s is the redshift space pair separation for the three different bias values given in Table[T] as well as b = 1 and b = 0.5 for comparison. 
The dashed lines show the predictions of the scale-dependent Gaussian streaming model, Eq. |25| which we evaluated using the linear theory expectations for 
fg( r )> vnCr), and <T^ 2 (r,[i). We see that accounting for the full real-to-redshift space mapping for a Gaussian pairwise velocity probability distribution that 
agrees with linear theory significantly modifies and £4 from the Kaiser formula expectation on quasilinear scales, with larger corrections at higher bias. 



Here r = r\ + y 2 is the pair separation in real space and ji = y/r. 
This model sums over all pairs with real space LOS separation y 
that are mapped to redshift space LOS separation r„ with a Gaus- 
sian probability distribution whose mean and dispersion depend on 
both scale r and angle with the LOS, \i. Note that our approach is 
not mathematically rigorous, but is rather an attempt to capture the 
most important non-perturbative effects. It can be thought of as a 
generalisation of models such as Eq.[23]to include the scale depen- 
dence of the exponential terms. 

Non-linear evolution causes the dark matter pairwise veloc- 
ity PDF to be significantly non-Gaussian on all scales; Eq. |25|ne- 
glects this. In Fig. [3] we compare the halo pairwise velocity PDF 
(solid) with that of the dark matter (dotted) for pairs separated by 
30 /r'Mpc < r < 31.5 /r'Mpc. The dashed lines show a Gaussian 
PDF with the same mean and variance as the halo PDF measured 
from the simulations. The PDF for the dark matter has a lower mean 
infall (as expected in linear theory). The dark matter PDF has 30 
per cent larger variance than the halos, and clear exponential tails 
at large velocities, as one would expect given that intrahalo virial 
velocities contribute only to the dark matter PDF. In contrast, the 
variance of the halo PDF is close to the expectation from linear 
theory. While skewness and kurtosis are clearly evident in the halo 
pairwise velocity separations at r = 30 /r'Mpc, if their scale de- 
pendence is smooth, £ ! may be relatively immune to them. There- 
fore, we will explore the accuracy of Eq.[25]for the redshift-space 
clustering of halos in the rest of this paper, and show that for many 
of the statistics of interest, this approximation is adequate on scales 
of tens of /r'Mpc. We also note that Tinker 1 2007} characterises 
the skewness and kurtosis of of the pairwise velocity PDF, and de- 
velops a physical model that describes them based on how a halo's 
environment depends on its mass. However, the improved accuracy 
of his model on smaller scales than are of interest in this paper 
comes at the cost of much greater model complexity than Eq. [25] 
His fig. 5 demonstrates that the relative deviations from Gaussianity 
quickly decline as a function of pair separation for r > 10 /r'Mpc, 
so we expect Eq.[25]to improve in accuracy as we move to larger 
scales. 

Allowing for the non-perturbative real-to-redshift space map- 
ping given in Eq. [25] alters of the redshift space correlation func- 
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Table 2. Constraints on /erg from Perc ival et al.| ^2004) , |Guzz o 
1 2007 1, da Angela et al. 12008 1, Samushia, Percival, & Raccanelli 
and |Blak^ret"aT"]p0*TTJ 7respectively. 



tion substantially in the quasilinear scales of interest in this paper 
(~ 30 - 80 /r'Mpc), even when we insert the linear theory pre- 
dictions for £ r g , v l2 (r), and cr^ /± (r) (Equations [7] to [IT]); we show 
this in Fig. [4] The m apping both smooths the BA~0 feature (but see 
also |Tian et al.|20lT| and enhances £2,4 on quasilinear scales, while 
only slightly suppressing Clearly, the amplitude of the correc- 
tion terms to the Kaiser prediction depend strongly on bias; we will 
study the bias dependence in more detail in Section[8] 



3 SETTING THE BAR: TARGET MODEL PRECISION 

Before entering a detailed study of the accuracy with which various 
models reproduce the clustering of halos in A'-body simulations, 
we consider the level of model accuracy desired in the near-term 
future. One aspect of this problem is the raw statistical power of 
current and upcoming surveys; for concreteness we focus on the on- 
going Baryon Oscillation Spectroscopic Survey (BOSS; Schlegel, 
|White, & Eisenstein|2009] |Eisenstein, et al.|2011^ . The other as- 
pect to consider is the level at which competing theories differ in 
redshift-space distortion observables, and on what scales the differ- 
ences are largest. 

For any expansion history H(z) and negligible dark energy 
clustering, general relativity makes a unique prediction for the 
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3.1 Fisher Matrix Analysis in Configuration Space 
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Figure 5. Fractional error on /cr 8 and berg as a function of the mini- 
mum configuration space separation .v„„„ used in the analysis, assuming 
a covariance matrix given by the usual Gaussian cosmic variance term and 
Poisson sampling variance term, with b = 2, V sur veY = 5h Gpc 3 , and 
n = 3.0 X 10~ 4 /! 3 Mpc~ 3 . We use the model in Eq. |23| and consider two 
cases. In the first, we marginalise over the value of cr, which controls the 
small scale damping. In the second case, we assume cr is perfectly known; 
in both the fiducial value of cr is 3.5 /r'Mpc. Marginalisation over cr in- 
creases the error on berg and /erg (upper two curves in both panels). The 
dashed curves show constraints when only fo,2 are used; the solid curves 
include £4 as well. In the case where cr is known (lower curves), C max = 2 
and 4 are indistinguishable. Finally, we compare the configuration space re- 
sults to an analysis in Fourier space (White, Song & Percival 2009), where 
a sharp cut-off in wavenumber is imposed, shown as dotted curves as a 
function of \.\5nlk max and also assuming cr is known. 



growth of structure as probed by redshift space distortions; alter- 
native theories of gravity and/or dark energy could potentially alter 
the predicted velocity field amplitude by ~ 10 per cent level (e.g. 
Song & Percival 2009 , but unfortunately there is no generic predic- 
tion for the level of deviation from ACDM+GR for alternative mod- 
els). Furthermore, combining precise measurements of peculiar ve- 
locities using redshift space distortions with lensing measurements 
provides an additional, powerful test of general relativity ('Zhang 
|et al.|2007|[Song & Koyama|2009| [Reyes et al.|2010| >. Therefore, 
extracting constraints from quasilinear scales where lensing can be 
measured with high signal-to-noise is desirable. 

In Table|2]we compile measurements of /cr 8 from galaxy red- 
shift surveys; see |Song & P ercival (20091 for the details on how 
these numbers were inferred from the reported constraints. In par- 
ticular, recent analyses of the SDSS LRG sample ( Samushi a, Per-| 
|cival, & Raccan e1li|201 l| l and from the WiggleZ Dark Energy Sur- 
vey ([Blake, et al. |201 \\ provide constraints on the growth of struc- 
ture at the ~ 7 per cent and ~ 5 per cent level, in the limit that 
each redshift bin can be treated independently. As we show in the 
next section, the Baryon Oscillation Spectroscopic Survey (BOSS; 
|Schlegel, White, & Eisenstein||2009[ |Eisenstein, eTaLl[2TJTT) will 
measure the peculiar velocity field at the few percent level. Thus, 
current and near-term observational constraints demand better pre- 
cision than is currently available from analytic models. 



In this section we present a calculation of the linear theory con- 
straints on berg, ferg, and cr, using Eq.[23]as our underlying model, 
and assuming perfect knowledge of P' m (k); the same calculation 
was done in White, Song & Percival ( 2009 1 in Fourier space. 

Throughout this section, we work at z = 0.55 where / = 0.74 
in our fiducial cosmology. Our fiducial value for cr is 300 km s -1 = 
3.5/r'Mpc. Fig. [5] compares the relative error on /cr 8 and berg 
as a function of the minimum redshift space pair separation s m j n 
used in the analysis; we assume s raax = 180/z~'Mpc. In the fig- 
ure, solid lines show constraints including £0,2,4, while dashed lines 
include only and £2- We consider two cases: one in which er 
is marginalised over (upper two curves in both panels), and the 
other where cr is assumed known (lower three curves in both pan- 
els). The parameter cr in Eq. |23| traditionally represents the effect 
of redshift errors and small scale velocity dispersions (fingers-of- 
god), but has also been used to absorb some of the theoretical un- 
certainty in the redshift space power spectrum (e.g., Percival & 
|White|2008[|Taruya, Nishimichi & Saito|2010) . If cr can be con- 
strained by better modelling of redshift space halo clustering in 
the quasi-linear regime and/or complementary observables sensi- 
tive the small-scale galaxy velocity distribution, Fig. [5] indicates a 
large sensitivity gain to the large scale peculiar velocity field am- 
plitude for s m j n ~ 30/i~'Mpc; in this paper we will focus only on 
the former approach. 

Except in the case of unknown cr, in the range 20/r'Mpc< 
■'mm < 40fr~'Mpc very nearly all of the information on /cr 8 is con- 
tained in £ and £>■ We therefore focus our modelling efforts on 
those two moments, though accurate modelling of £ 4 will be impor- 
tant if the geometric factors At and H are also being fitted 1 Taruya,| 



Saito & Nishimichi 2010; Kazin, Sanchez & Blanton|201 1^ 7 

Finally, in the case of fixed, perfectly known cr, we can com- 
pare the constraints in configuration space to those in Fourier space 
using the public code released with White, Song & Percival ( 2009 1, 
shown as dotted curves as a function of 1. 15 n/k max , which provides 
an empirical mapping between s mm and fe max that gives nearly iden- 
tical errors on /cr 8 over the range of scales we consider. Reassur- 
ingly, the same fractional errors are recovered in both configuration 
and Fourier space when all scales are included. While the relatively 
smaller errors on berg in Fourier space indicate that there is more 
information on the shape of the power spectrum at fixed non-linear 
scale .s m j„ ~ n/k max , the present analysis ignores the fact that shot 
noise will contribute at all k to the power spectrum, and cannot be 
assumed Poissonian at this level of accuracy for halos or galaxies 
(e.g. |Cooray & Sheth|[2002l |Smith, Scocciimarro, & Sheth||2007l 
|Hamaus et al.|2010| >. 



4 COMPARING THE SCALE-DEPENDENT GAUSSIAN 
STREAMING MODEL FOR HALOS WITH N-BODY 
SIMULATIONS 

4.1 Simulations 

To assess the accuracy of our analytic model we compare it to 
halo correlation functions derived from the 20 Af-body simulations 

1.5 /j^'Gpc and parti- 
~ for 



oflWhite et al.U201l]) which have L b , 



White et al. 



2011 



cle masses m p = 7.6 x I0 >0 h [ M Q (see 
more details). The simulations are of a ACDM cosmology with 
= 0.274, n A = 0.726, h = 0.7, n = 0.95, and cr 8 = 0.8. In this 
section we consider only a single snapshot at z = 0.55. The halo 
catalogues are generated by the Friends of Friends (FoF) algorithm 
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Figure 6. Redshift space Legendre moments fo,2,4 as a function of red- 
shift space separation s divided by the linear theory expectations. Points 
with eiTors show measurements from A'-body simulations for the high (cir- 
cles), low (triangles), and HOD (diamonds) samples described in Table [T] 
In the upper left panel we offset the "high" and "low" halo sample by ±10 
per cent, respectively, for clarity. Black solid curves show the result of the 
scale-dependent streaming model, Eq.[25] in the case where fjj(r), vi 2 (r), 
and <r 2 l2 (r,fi) have been measured directly from the A'-body simulation real 
space clustering and velocity statistics for the same halos. This comparison 
tests the approximations in the form of Eq. |25| supposing we had a perfect 
theory for the input real space statistics. The bottom right panel shows the 
fractional error of the streaming model on £2 as a function of scale. For all 
halo mass bins we studied, the streaming model is accurate at the 2 per cent 
level (or better) at s = 25 /r'Mpc. 

( |Davis et al.|1985| > with linking length 0.168 times the mean inter- 
particle spacing. The principal purpose of this investigation is to 
generate an accurate model of redshift space distortions as a func- 
tion of halo bias. Throughout this section we will compare models 
with three distinct halo samples; the same level of agreement was 
found for all mass bins we investigated. The detailed properties of 
the samples are listed in Table [T]- the 'high' and 'low' samples 
impose sharp halo mass cuts, while the 'HOD' sample includes 
all halos hosting a mock galaxy in the catalogues used in |White| 
|et al.| |201 1 ). This catalogue is different from the mock galaxy cat- 
alogue only in that the "satellite" galaxies are excluded from the 
sample. The purpose of this sample is to demonstrate that a sample 
including a broad range of halo masses can still be described by our 
model, bringing us one step closer to modelling what we observe: 
the galaxy density field. 

4.2 Results 

In this section we test the approximations entering the scale- 
dependent Gaussian streaming model, Eq.[25] for halos in A'-body 
simulations. We do so by assuming we have a perfect model for 
all the real-space quantities entering Eq. [25] the halo correlation 
function £' h (r), the mean pairwise infall velocity as a function of 
separation, v 12 (r), and the pairwise velocity dispersions parallel and 
perpendicular to the pair separation vector, crj ± (r). We simply mea- 



sure these quantities from our A'-body simulations and input them 
into Eq. [25] to generate model curves for Fig. [6] shown as black 
solid curves (we will consider analytic models for these functions 
in Section [5~T| l. 

Fig. [6] shows £0,2,4 f° r the high (circles), low (triangles), and 
HOD (diamonds) halo subsamples in TableJT] these are the same as 
in Fig. [2] We again offset the 'high' and 'low' samples by ±10 
per cent for clarity in the £ panel. In the lower right we plot 
the fractional difference between the simulation results and our 
model for the quadrupole £ 2 , highlighting the percent level agree- 
ment between the streaming model and the A'-body results down to 
s ~ 25 /r'Mpc. 

For clarity we omit the results on BAO scales; the agreement 
is as good as the Matsubara (2008b I model shown in the bottom 
right panel of Fig. [2] On smaller, "quasi-linear" scales, we find a 
drastic improvement over the model of Matsubara ( 2008b I, which 
is unsurprising since our model includes real space halo cluster- 
ing and velocity statistics derived from A'-body simulations, and 
treats the non-linear mapping between real and redshift space non- 
perturbatively. The streaming model provides a sufficiently accu- 
rate description of the monopole <fo down to 10 /r'Mpc, with no 
detectable deviation even for this large volume of simulations. The 
streaming model is an accurate model for f 2 at the 1-2 per cent level 
for redshift space separations greater than ~ 25 /r'Mpc. However, 
both the LPT and streaming model predictions are poorer descrip- 
tions to the A'-body results for higher multipole moments. On scales 
s < 80 /r'Mpc, £4 deviates from the linear theory prediction by 
0(1). However, the streaming model has the right dependence on 
bias and may be sufficiently accurate on some range of scales for 
upcoming surveys, given the large measurement errors on £4. We 
are also optimistic that treatment of non-Gaussian small-scale ve- 
locities that are uncorrelated with the linear theory pairwise infall 
velocities we wish to measure may absorb some of the remaining 
discrepancy. We reserve this line of study for future work. 

As we emphasised earlier, both P m and P m are suppressed 
compared with linear theory on relatively large scales (k > 
0.03 &Mpc~'), and we expect that suppression to affect redshift 
space clustering on similar scales. Even at 5 = 50 /r'Mpc, £ 2 is 
suppressed by 2, 6, and 8 per cent compared to the linear the- 
ory expectation, in line with the findings of Oku mura & Jing| 
< |20 1 1| >. Comparison with Fig.[4]shows that on the scales of interest 
(~ 30- 80 /r'Mpc), both non-linear corrections to the velocity field 
and the non-linear mapping from real to redshift space are relevant 
to predicting £ 2 , at a level comparable to the latest ~ 7 per cent pe- 
culiar velocity field constraints ( [Samushia, Percival, & Raccanelli| 
|201 1| |Blake, et al.|20li] >, and certainly at the level of constraints 
projected for BOSS in Fig.[5] 



5 COMPONENTS OF THE SCALE-DEPENDENT 

GAUSSIAN STREAMING MODEL: PERTURBATION 
THEORY 

In this section, we consider perturbation theory models for the com- 
ponents of the scale-dependent Gaussian streaming model in pur- 
suit of an analytic model for the redshift space correlation function 
of biased tracers. While our model for the real space halo correla- 
tion function ^' h {f) includes non-linear Lagrangian bias, our pertur- 
bation theory calculations for halo infall velocities and dispersions 
were carried out in standard perturbation theory, where consider- 
ing non-linear biasing requires the definition of a smoothing scale 
( |Heavens, Matarrese, & Verde|1998"l >. |Roth & Porciani| l |201 \) also 
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5.2 Mean halo infall velocities, V[ 2 (r) 




A mean (pairwise halo) velocity along the pair separation vector 
arises from the correlation of the density field with the velocity 
field: 

<[1 + b5(x)][l + b6(x + r)][v(x + r) - v(x)]> 



v l2 (r)r ■■ 
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Figure 7. Real space correlation functions for the halo samples in Table 
[T] The high (low) bias bins are offset by ±10 per cent for clarity. The solid 
curve shows the prediction of Matsubara 1 2008b) when we fit for the large 
scale bias as a free parameter. For the dashed-dot curves, we show the Mat- 
|subara| (2008b) prediction when we artificially set b 2 = 0, just to show the 
order of magnitude of the contribution from non-linear halo bias. Using the 
linear theory covariance matrix, we find a good fit for scales r > 25 /i~'Mpc 
to our simulation results, which total 67.5 /i~ 3 Gpc 3 . BAO scales are not 
shown here, but the fit is good on those scales as well. 



found that the Eulerian local biasing scheme is not very accurate, 
and it is not equivalent to a Lagrangian biasing scheme ( Matsubara] 
|201 1| >. For these reasons, we perform our velocity calculations for 
the simplest model of linearly biased tracers with bias b, which we 
identify with the large scale bias we fit to the real space correlation 
function using LPT. 



5.1 Halo real space correlation function, f£(r) 

The Lagrangian perturbation theory prescription of Ma tsubaraj 
(2008b) for describing local Lagrangian biased objects includes 
both first and second order bias terms, which are related by the 
peak-background split to the halo mass function. This theory pro- 
vides a good description (accurate at the 1 per cent level) on 
scales r > 25/i~'Mpc. We show this explicitly for scales r < 
80/r'Mpc in Fig. [7] and note that the fit is good on the full range 
of scales we study (r < 180/z~'Mpc). If we fit the LPT predic- 
tion to £; r h {r) for separations 30/i~'Mpc < r < 180/i~'Mpc, we 
find x 1 = 93, 104, 119 for 99 degrees of freedom and one free 
parameter (the large-scale halo bias) for the high, low, and HOD 
halo subsamples, respectively. Fitting instead to linear theory gives 
X 1 = 212, 274, 403. In both cases we use the standard linear theory 
with the Poisson sampling assumption to derive the covariance ma- 
trix (Af(rf)Af(r/)V To illustrate the amplitude of the second order 
bias corrections for the halos of interest, we also plot the LPT pre- 
diction when bi is artificially set to zero (dashed-dot curve) com- 
pared with the LPT prediction including nonzero b 2 (solid curve); 
the second order bias contribution to ff is quite small for the linear 
halo bias values we consider. 



<[l + M(x)][l+fcrJ(x + r)]> 



(26) 



where b is the linear halo bias. In perturbation theory, the density 
and velocity fields are written as a sum of terms (6 = S\ + 62 + 
Sj, + ...), with the subscript denoting the order of their dependence 
on the linear density field, 5i(k). Up to fourth order in 6i(k), there 
are three distinct corrections to the linear theory expectation V\ 2 {r) 
given in Eq.|7] each with a different dependence on bias: 

[l + b 2 ? m {r)\ v^(r)f = 2Z><<5 1 (x)v 1 (x + r)> + 

2b ^ <<5,(x)v4-,(x + r)> + 2b 2 J] (<S,(x)<S,(x + r)v 4 -, w (x + r)) . (27) 

i>0 i,j>0 

The three (<5,v<4_,) terms arise from the perturbation theory correc- 
tions to P 6e , and the bias dependence is the same as the linear theory 
term, (S^'i). The terms from three-point correlations (S{SjV^-i~n 
scale with b 2 , so their contribution will be larger for more highly 
biased tracers. Note that these terms are exactly the ones evaluated 
in Appendix B of |Tang, Kayo, & Tak ada (2011). We provide ex- 
plicit expressions for all of these terms in Appendix [A] Finally, the 
pair- weighting correction, 1/[1 + b 2 ^ n (r)], will be larger at a given 
scale for more biased objects. The relative contribution for these 
three corrections is shown in Fig. [8] for b = 2. At least for b = 2 
of interest to BOSS, the two-point corrections from P se never dom- 
inates, so only including the two-point corrections (as in Eq. [24| 
will be a poor model for the redshift space power spectrum; we 
should expect important contributions from the bispectrum as well 
(Scoccimarro 2004; Taruya, Nishimichi & Saito 2010; Tang, Kayo, 
|&Takada|2011| >. 

In the left panel of Fig. [5] we compare the deviations from 
linear theory infall velocity predictions measured from our simu- 
lations to our perturbation theory calculation. The expected V\ 2 (r) 
depends on the halo bias, and for this we use the first order bias de- 
duced from fitting the real-space halo clustering to the LPT model 
of IMatsubara] p0~0 8b) ; Table [T] indicates that the best fit LPT bias 
can differ at the few percent level from the best fit linear bias. At 
the percent level, Fig.|9]shows that the LPT bias predicts the correct 
infall velocity amplitude on the largest scales. This confirms the 
common assumption in the literature that "velocity bias" is small, 
at least for halos in the bias range we have studied. Perturbation 
theory provides a relatively good description of the departure from 
linear theory. The difference between the simulations and perturba- 
tion theory depend on halo bias and agree best for the HOD halo 
subsample, which Fig. [7] indicates is the sample with the smallest 
second-order bias. We note that there is good theoretical motivation 
to expect the bias relevant to the matter-velocity cross-correlation 
to differ from the one inferred from clustering, and have scale- 
dependence (Desjacques & Sheth 2010). While including a more 
complicated biasing model may improve agreement with the pair- 
wise velocity statistics of halos, our results suggest that these cor- 
rections are small on the scales of interest here. 



5.3 Halo velocity dispersions, cr 2 2 (r,i/) 

Analogous to Eq. [27] we compute the pair-weighted velocity dis- 
persion. The result depends only on fi 2 = cos 2 (0 fr ), where <p (r is 
the angle between the LOS and the pair separation vector. As in 
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Figure 9. Left panel: Pairwise halo infall velocities for the halo subsamples listed in Table [TJ divided by the linear theory expectation (Eq.|7J; the high 
(circles) and low (triangles) bias bins are offset by ±10 per cent for clarity relative to the HOD sample (diamonds). Note that the expected vj?(r) depends 
on halo bias; we use the best fit LPT bias given in Table [JJ The solid curves show the prediction for linearly-biased tracers from our perturbation theory 
calculation, Eq.[^] Middle and right panels: Pairwise halo velocity dispersions measured from the simulations are shown as dotted curves for the directions 
along and perpendicular to the pair separation vector. For cr 2 (r) we subtract the square of the mean infall Vj 2 (r) to get the dispersion about the mean, and 
account for this accordingly in the perturbation theory prediction. On large scales there is a small offset from the linear theory predictions (dashed curves). We 
shift the measured dispersions by a constant value (lower dotted curves) to compare the scale-dependence with both linear theory (dashed) and perturbation 
theory (solid; also shifted to agree with linear theory at r = 180/i~'Mpc). For <x 2 (r), the agreement between the perturbation theory calculation and the halo 
dispersions in simulations is so good that the curves are difficult to distinguish in the figure. 



Eq.[TT| we present results for the velocity dispersion perpendicular 
and parallel to the LOS, which can be combined to give the disper- 
sion as a function of (r,fi 2 ). 



((1 + b8(x))(\ + bS(x + r))(v ( (x + r) - v'(x)) 2 
<(1 + bS(x))(\ + bS(x + r))) 



(28) 



We again separate the terms by their bias and scale dependence, 
and provide explicit expressions in Appendix|A| 

[l+b 1 ar)]<T] 1 (r, l i 2 ) = 

2(((/(x)) 2 - Ax)v f (x + r)))+ (29) 

2b(S{x)(v'(x)) 2 ) + (30) 

2b (((5(x)(/(x + r)) 2 } - 2 (<5(x)/(x)/(x + r)}) + (31) 

+2b 2 [(S(x)6(x + r)(v f (x)) 2 } - (S(x)S(x + r)v f (x)v f (x + r)}] . (32) 

The higher order terms in line|29]can be accounted for in the form 
of Equations |9] and [10] but replacing P' m {k) with the perturbation 
theory result for Pgg{k). The term in line [30] evaluates to a con- 
stant, which for our fiducial cosmological parameters is 13.24/?/ 2 
(/r'Mpc) 2 . While the perturbation theory calculation overestimates 
the amplitude of the effect, Fig. [9] does show a large scale off- 
set between the linear theory velocity dispersions and those mea- 
sured for halos in our simulations; the offset for our HOD sample 
is ~ 1.5 (/i~'Mpc) 2 . There is a slight dependence on halo mass, 
with the high mass sample offset ~ 1.0 (/j~'Mpc) 2 and the low 
mass sample offset ~ 2.2 (/r'Mpc) 2 . The offsets for the parallel 
and perpendicular components are in agreement with each other 
at the level of 0.1 (/r'Mpc) 2 . In a future paper we expect to ac- 
commodate this offset into our theory along with other small-scale 
isotropic dispersions due to redshift errors and fingers-of-god from 
satellite galaxies. However, we point out that for halos, the small- 
scale, incoherent velocity dispersion is much smaller than the linear 
theory dispersion cr 2 = 21 (/r'Mpc) 2 that appears in Eq.|24| 

Line[3T]has many terms that depend on B me \ we evaluate them 



explicitly in the Appendix[A] Finally, line[32]reduces to 

b^Jr^-Jr^ 1 ) + \v\ 2M {r)n 2 (33) 

This shows that the pair-weighting factor cancels out to leading or- 
der, and the dispersion should be increased along the separation 
vector due to the linear infall. The middle and right panels of Fig. [5] 
compare the scale dependence of the HOD halo subsample veloc- 
ity dispersions parallel and perpendicular to the LOS, respectively. 
The upper dotted curves show the dispersions about the mean infall 
measured from the halos in the simulations. To compare the scale 
dependence, we subtract a constant from the measured velocity dis- 
persions (lower dotted curves), forcing agreement with the linear 
theory expectation (dashed curves) on the largest scales. The solid 
curves show the perturbation theory expectation after subtracting 
the expected mean infall contribution, {vn(r)) 2 n 2 , and a constant 
to force agreement on the largest scales. As the figure shows, the 
agreement between the scale-dependence of the dispersions pre- 
dicted from perturbation theory and the simulations is excellent for 
the HOD subsample. We find that it is only slightly worse in the 
radial direction for the other mass bins. 



6 ACCURACY OF THE PERTURBATION THEORY 
SCALE-DEPENDENT GAUSSIAN STREAMING 
MODEL 

By combining the results in Section [5~T] using Eq.[25] we have an 
analytic prediction for £0,2,4- Because £ 2 contains almost all of the 
available information on /cr 8 (see Fig. [5}, we focus on that quan- 
tity here. Fig. [10] shows the same simulation results as in Fig. [6] 
for £ 2 (symbols with errors) while the lines show the prediction of 
the perturbation theory scale-dependent Gaussian streaming model. 
While the agreement is excellent for the HOD halo subsample, the 
accuracy of the model is worse than 2 per cent for scales smaller 
than ~ 40 h Mpc. This disagreement can be traced almost entirely 
to insufficient accuracy in the perturbation theory prediction for 
V\2(r). Replacing the perturbation theory prediction for v 12 (r) with 
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Figure 8. The three components of the next-to-leading order correction 
to the linear pairwise infall, v^"(r); each term scales differently with bias. 
The solid curve shows the total prediction for b = 2 halos at z = 0. 55 
The dashed curve shows the contribution from the (SjVj^ terms in Eq. 27 
(i.e., those originating from corrections to Pse), the dotted curve shows the 
contribution from b (SiSjVkj terms (i.e., those originating from Bsss), and 
the dash-dot curve shows the pair- weighting correction contribution, 1 / [l + 
b 2 ^ r m (r)] - 1. Note that the Pse correction terms never dominate the total 
correction, so we should generically expect that models neglecting the three 
point and pair- weighting corrections to be poor models of the redshift space 
power spectrum for halos (e.g., Eq.|24). 



the simulation results, keeping all other terms fixed, dramatically 
improves the agreement with simulations. This explains why £ 2 for 
the HOD model agrees so well for the HOD subsample - it is for 
this sample that the perturbation theory model for v 12 agrees best 
with the simulations. Eq. [T7] demonstrates that even in the linear 
limit, both v l2 (r) and its derivative contribute to f 2! s ° the degra- 
dation in accuracy from vi 2 (r) to £ 2 is primarily because the slope 
of vJJ(r) is inaccurate for the low and high mass halo sample, but 
coincidentally correct for the HOD subsample. We plan to explore 
other perturbation theory approaches for v I2 (r) in a future paper in 
hopes of improving the accuracy of the model. 
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Figure 10. Inserting the results of Section [5~71 into Eq. |25| gives an analytic 
model for f 0.2,4 ■ We compare the f 2 prediction of this model (solid curves, 
with lower bias predicting lower f 2 /f 2- ii n ) to the same simulation results as 
in Fig. [6] The fit is excellent for the HOD subsample. The inaccuracy for the 
other mass bins can be traced to the insufficient accuracy of the perturbation 
theory prediction for vi 2 (r) and its slope on small scales. 



example to emphasise that the non-linear corrections are not nec- 
essarily smaller at higher redshift, since more highly biased objects 
are often being selected. We have checked the behaviour of the per- 
turbation theory predictions at z = and z = 1 for the real space 
statistics examined in Section |5~T) and find good agreement with 
naive expectations. At fixed b = 2.8, the sudden upturn in the LPT 
prediction for t; occurs at increasingly large r as structure evolves 
(22, 26, and 31 /r'Mpc at z=l, 0.5, 0, respectively). At fixed z, the 
upturn increases with b, as can be seen in Fig. [7] We find similar re- 
sults when comparing pairwise halo velocity statistics; in general, 
the perturbation theory predictions are better at higher redshift and 
lower bias. However, the "sweet spot" in the perturbation theory 
prediction at b = 1.9 in Fig.|9]persists at z = as well, though the 
halo mass range at fixed b depends on z- 



7 REDSHIFT DEPENDENCE 

Though we have focused on z = 0.55, Fig.[TT]shows that the scale- 
dependent Gaussian streaming model works extremely well for 
massive halos both at z = 1 and z = 0. In the upper panel, we com- 
pare samples containing all halos above M min = 10 13387 /i~'Mq, 
for which the best fit bias values are bipr = 3.83,2.79 and 1.88 
(in order of decreasing redshift). In the lower panel, we vary 
M mill = 10 129 , 10 13 - 387 , 10 1395 /j-'Mq, so that the bias of the sam- 
ples is fixed at b LPT ^ 2.8. The redshift dependence of £ 2 /f 2 .ii n 
on small scales is opposite in the two panels, which can be under- 
stood simply in terms of the absolute real space clustering ampli- 
tude that enters the convolution in Eq. [25] At the fixed value of 
M we have chosen, f r oc b{z)cr i {z) decreases as structure grows, 
whereas at fixed bias, f r increases with time as <x g (z). We use this 



8 BIAS DEPENDENCE 

As we have emphasised, there are two relatively large corrections 
to the Kaiser prediction for <f 2 . In this section we illustrate explic- 
itly how they depend on bias. The suppression of the halo velocity 
statistics relative to linear theory lowers the amplitude of £ 2 , while 
the non-linear mapping between real and redshift space increases 
the amplitude of £ 2 in the bias range we have studied (see Fig. [4). 
To illustrate the effect of the former, we use the Kaiser limit map- 
ping between real space clustering and velocity statistics (Equa- 
tions [16] through ITS), but input the real space W-body simulation 
results for these quantities rather than the linear theory expecta- 
tions. The results for the halo subsamples in Table[TJare the lower 
three curves in Fig. [12] The weak dependence of the non-linear 
corrections on the velocity statistics (Fig. [9) translates into a rela- 
tively weak dependence on halo bias on £ 2 . In contrast, the effect of 
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Figure 11. To check the redshift dependence of our model, we show 
%zl%1,lin at different simulation snapshots z = 1 (circles), z = 0.55 (tri- 
angles), and z = (diamonds). In the upper panel, we select all halos above 
a fixed minimum halo mass M„„„ = 10 13 387 /t'Mq. In the lower panel, we 
vary the minimum halo mass included in the sample such that the sample 
bias remains fixed at 2.8; the z = 0.55 sample (the 'high subsample' in pre- 
vious plots) is the same in both panels. The solid lines are the predictions 
from the Gaussian scale-dependent streaming model with input real space 
clustering and velocity statistics measured from the W-body simulations (as 
in Fig.[6). 



the non-linear mapping depends strongly on halo bias. To see this 
explicitly, one can expand Eq.[25]by assuming that the pairwise ve- 
locity PDF l P(v : ;r) is a smooth and slowly varying function of r; 
this procedure will be more accurate at small /i, where a smaller 
range of real space separations contribute pairs at a given redshift 
space separation. Eq. 53 of Scoccimarro (2004) does this in the case 
of the exact Gaussian result, and the same terms (along with many 
others) appear when the expansion is performed on our Eq.[25] We 
have verified that the dominant non-linear correction term for £0,2 
in our bias range comes from the term -d / dy[t;v l2 ]: 

Nimappmgin = (H ~ U " /" ~T ■ ( 34 ) 

The upper curves in Fig.[T2]are the same as the lower ones, but with 
this extra term included to approximate the non-linear mapping 
step; these predictions are in reasonable agreement with Eq. [25] 
but performing the full integral is a noticeably better fit to the sim- 
ulation results. What we wish to emphasise is that the non-linear 
mapping produces a term (Eq. [34} that contributes to £ an d f 2 
and scales like b 3 . This is in disagreement with the recent results 
of |Tang, Kayo, & Takada| ( |2011) , who use a non-linear correction 

We 



27 



term equivalent to the fijj6fV4-j- A contribution in our Eq. 
can see why that term (the dotted curve in FigureJH} provideTa rea- 
sonable fit to their simulation results at one value of b: its shape 
roughly mimics our non-linear mapping term that dominates on 
small scales. However, our more detailed analysis demonstrates 
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Figure 12. The predictions for £ 2 for halos (b=2.8, solid; b=1.9, dotted; 
b=1.4, dashed) for two different mappings between real and redshift space 
using the simulation results for the non-linear real space clustering and ve- 
locity statistics. The first mapping (lower curves at s = 30/i _1 Mpc) is equiv- 
alent the mapping assumed in the Kaiser formula (our Equations | 1 6| through 
|18( , while the second includes the dominant correction term (Eq. |34) from 
the non-linear mapping, which scales like b 3 . 



that many other terms are of comparable size to the one considered 
in |Tang, Kayo, & Takada|p0TT) . 

|Okumura & Jing (2011) find that the value of yS recovered 
from massive halos b > 1.5 is relatively close to the expected lin- 
ear value, but lower mass halos recover a smaller value compared 
with linear theory. Fig. |12| illustrates why: for our central galaxy 
sample (dotted curves), the non-linear effects of velocity suppres- 
sion and real-to-redshift space mapping approximately cancel for 
s > 30/r'Mpc; at low halo bias, we expect the non-linear mapping 
corrections to be small, and the measured £ 2 should be closer to the 
lower curves. Of course, the bias where this near-cancellation oc- 
curs will depend on redshift, and because of the Z? 3 correction term, 
it will only be true in a limited range of bias values. 

To be more quantitative, for the halo bias range we have stud- 
ied (b = 1 .4 - 2.8), fitting the Kaiser formula to <f an d £2 to derive 
constraints on b and / on scales r ~ 30/r'Mpc will bias the con- 
straints on / by +2, -6, and -10 per cent for b\ m = 2.67, 1.84, 1.41, 
respectively, under the assumption that the smallest scales included 
dominate the signal-to-noise. These biases are already at the level 
of current statistical errors. 

While our simulations cannot reach halos with b ~ 1, our 
analysis can shed some light on what behaviour to expect for 
b ~ 0.8 - 1.2 halos of relevance to the WiggleZ survey {Blake, et] 
|al.|2011| > and more closely related to the perturbation theory stud- 
ies for matter. Fig. [4] shows that the non-linear mapping should be 
a small correction for r > 30fr~'Mpc, but that it still amplifies £ 2 
on smaller scales. If we evaluate our perturbation theory predic- 
tions for b = 1, we find that the total correction to v 12 (r) can be 
well-approximated by only the P 6tj term; however, the higher or- 
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der terms still dominate the corrections to the velocity dispersions. 
Therefore, while the bispectrum terms should not be negligible, fit- 

with a 2 v treated as a free parameter 



24 



ting formulae based on Eq. 
capture at least some of the relevant non-linear corrections and can 
absorb the rest into cr; ; even while providing a good fit to the data, 
we are not guaranteed that the underlying peculiar velocity field 
amplitude will be recovered at the few per cent level, especially 
when fitting down to k max = 0.3 h Mpc~' (i.e., s„„„ » 10/r'Mpc). 



9 DISCUSSION 

In contrast to many recent theoretical investigations of redshift 
space distortions, which have focused on the matter density field 
and/or performed analyses in Fourier space, in this paper we focus 
on the redshift space clustering of dark matter halos in configura- 
tion space, and use 67.5 /T 3 Gpc 3 of ./V-body simulations to make 
precise measurements of £ ,2,4 as a function of halo bias. In our 
modelling we focus on two distinct corrections to the linear the- 
ory predictions: the non-linear mapping between real and redshift 
space, and the non-linearity of the halo pairwise velocities. We find 
both corrections to be important on the quasilinear scales of interest 
to this work (~ 30 - 80/r'Mpc). 

To model the non-linear real to redshift space mapping of 
pairs, we take a non-perturbative approach and employ the scale- 
dependent Gaussian streaming model, Eq. [25] in which the pair- 
wise velocity probability distribution function (PDF) is assumed to 
be Gaussian, but where the pairwise velocities have a mean and dis- 
persion that depend on the pair separation distance r and the angle 
of the pair separation vector with the LOS. A similar model has 
been used with some success on somewhat smaller scales ( |Tinker,| 
|Weinberg, & Zheng|2006l |Tinker|20"07) ) in the context of the halo 
model. Fig.[4]shows that this model significantly enhances both £ 2 ,4 
on quasilinear scales for real space statistics expected in linear the- 
ory. When we know perfectly the real space clustering and velocity 
statistics, Fig. [6] shows that this model is accurate at the < 2 per 
cent level for s > 25 /r'Mpc, i.e., at about the level demanded by 
ongoing experiments like BOSS. 

For the first time (to our knowledge), we have computed the 
next-to-leading order corrections to pairwise mean infall velocities 
and dispersions for linearly biased halos as a function of real space 
separation in standard perturbation theory. While we find relatively 
good agreement between our calculations and halo samples drawn 
from a large volume of N-body simulations (see Fig.|9j, when used 
as input into the streaming model, there remain offsets at the sev- 
eral percent level for s < 40fr~'Mpc. We are able to trace the 
discrepancy to the scale-dependent inaccuracy of the perturbation 
theory prediction for halo infall velocities; note that even in linear 
theory, the redshift space quadrupole £ 2 depends on the derivative 
dvn(r)/dr. In future work we hope to explore other perturbation 
theory schemes as a means to improve the level of accuracy of the 
mean pairwise velocity perturbation theory prediction; more com- 
plex biasing schemes may also improve agreement (Desjacques & 
|Sheth|2010| i. 

Both of the corrections we described above have significant 
higher order contributions that depend on the halo bias. In Section 
[3] we show that while the bias dependence is detectable but weak 
for the non-linear velocity corrections, the dominant correction for 
the non-linear mapping from real and redshift space scales as b 1 . 
These effects have opposite signs in the halo bias range we have 
studied, so that for some limited range of halo bias and scale, they 
approximately cancel. Our findings demonstrate, in line with other 



recent works, that a model of the form in Eq.[24]fhat only includes 
two-point corrections cannot accurately describe the dependence of 
redshift space halo clustering on bias. 

Finally, we note that in order to infer /cr 8 from redshift space 
distortions in halo clustering, one must make the assumption that 
the bias inferred from real space clustering is the same one that 
determines the halo pairwise infall velocity amplitude. Our large 
volume of ,/V-body simulations allows us to confirm this assumption 
at the per cent level on large scales, once we incorporate the effects 
of non-linear growth using perturbation theory. 

Before our model can be applied to analysing real galaxy sur- 
veys, the large velocity dispersions of satellite galaxies must be 
accounted for separately in the model. We are cautiously optimistic 
that these additional corrections, at least for will be relatively 
small. In a preliminary study, we found that £; 2 ( s = 30/r'Mpc) is 
damped by < 2 per cent when satellites from the best fit HOD of 
|White et aT] ( |201 1) are included. However, our model prediction for 
£4 is worse than for the halo-only samples. We reserve this line of 
research for a future work. 
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APPENDIX A: PERTURBATION THEORY 
CALCULATION DETAILS 



We refer the reader to Bernardeau et al. ( 2002| for an introduction 
to standard perturbation theory. In this appendix we provide explicit 
expressions for all terms contributing up to fourth order in the linear 
density field <5i(k) to pairwise mean velocities and dispersions for 
linearly biased tracers. Our results are presented using the Fourier 
convention given in Eq.|12| 



Al Mean halo infall velocities 

Terms in Eq.[26]that depend on only one 6 can be written in terms 
of the perturbation theory density-velocity cross power P sg : 



2b 



<<5i(x)vi(x + r)> + <<5 i (x)v 4 _,(x + r)> 

;>o 



dkkP%(k)Mkr) 



(Al) 



We evaluate P?J using the publicly available Copter code 1 Carlson, 



|White, &T admanabhan 2009 1. Strictly speaking, the perturbation 
theory result for Pse(k) becomes negative at large k; we simply trun- 
cate the integral over k at k max = 2 h Mpc~', where linear theory 
predictions are already well converged. 
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There are three distinct contributions of the form V4_j_y): 
<<Si(x)<S 2 (x + r)vi(x + r)) = 
2 I — ^P(|k 2 +k 3 |)/'(k3)F 2 (-k3,k 2 +k3y^ +k ^' r (A2) 

<<5i(x)<Si(x + r)v 2 (x + r)) = 



; rd 3 k 2 d 3 k 3 ij 
'J (27r)«fc 2 

, r d 3 k 2 d%ij 

J (2nfk\ 



P(|k 2 + k 3 |)P(k 2 )G 2 (-k 2 ,k 2 + k 3 y 



(k 2 +k 3 )T 



(S 2 (x)S l (x + r)v t (x + r)} 



, rd%d%ifk 3 

J (2nfk\ 



/>(k 2 )/>(k3)F2(k2,k3K 



(k 2 +k 3 )T 



(A3) 



(A4) 



where F 2 and G 2 are the standard perturbation theory kernels for 6 2 
and # 2 , respectively. 

Performing some of the angular integrations, we find 

<(5i(x)(5 2 (x + r)vi(x + r)> + (6i(x)6 { {x + r)v 2 (x + r)> = 



f*_ 
4tt 4 



Jo Jo 



du dy u P(u)P(yu)ji(ur) X 



2y(-19 + 24y 2 - 9y 4 ) + 9(-l + y 2 ?T 
168y 

(<5 2 (x)(5i(x + r)vi(x + r)> = ^ j x (ur)du x 

3yw - 10yw 3 + lw 2 



P(uy)dy J dwP(u^[l 



+ y 1 — 2yw) 



14(1 +y 2 -2yw) 



(A5) 



(A6) 



where we have defined T = tanfr'(2y/[l + y 2 ]). Note the last 
integrand can also be expressed as a sum of products of one- 
dimensional integrals that each depend on r. We prefer the ex- 
pression above, which can be combined with Eq. |A5| before do- 
ing the final integral over j\(ur). We caution the reader that 
the integrals are only well-behaved in combination; for instance, 
<<5i(x)5 2 (x + r)vi(x + r)> and (6\{x)8\(x + r)v 2 (x + r)> do not con- 
verge individually for a CDM-like power spectrum, where P(k) ~ 
k~ 2b at large k. Moreover, the integral over x = uy in Eq. A5 
yields a constant at large u, which shou ld then be integrated against 
u 3 P(u)ji(ur); combining with Eq. A6 before performing the inte- 
gral on u yields in integrand that behaves like u~ 9 j t (ur) at large u, 
which converges. 



A2 Halo velocity dispersions 

In this section we evaluate the contributions to <x 2 2 (r,/j 2 ), given in 
Equations[28]to[32] We first introduce some shorthand notation: 



J(fi (r , kr) = n (r Jo (kr) + (1 - |/ fr ) - 



kr 



kr 



h,m(y) = f siri"(e)cos"\G)e iycom dG. 
Jo 



(A7) 
(A8) 



The higher-order contributions to line|29]can be expressed in terms 
of the perturbation theory result of the velocity divergence power 
spectrum: 



(r (x)v'(x + r>) = ^ J dkP%(k)JQii,kr), (A9) 



1L 

In 2 



which is equivalent to Equations[8]-[To]at linear order. We evaluate 
Pm with Copter as well, again truncating the integral over k at 
k max =2h Mpc~\ where linear theory predictions are already well 



converged. The term in line 30 adds a constant to <r 2 2 (r, fi 2 r ) 



(<5(x)(v f (x)) 2 } = 2(<5,(x)^(x)v 2 (x)) + (<S 2 (x)(^(x)) 2 } 



6w 4 



s dk dy P(k)P(ky) x 



2y(3 + 8y 2 - 3y 4 ) + 3(-l + y 2 ) 3 T 
56y ' 



(A10) 



where again neither term converged individually. There are five dis- 
tinct contributions to the terms in line|3 1 1 

((<5(x)(/(x + r)) 2 } - 2 (<5(x)/(x)/(x + r))) = 

2(6 l (x)v' l (x + r)v e 2 (x + r)) + (<5 2 (x)(v^(x + r)) 2 } 

-2 (^,(x)v^(x)v^(x + r)) - 2 (^(x^Cx^Cx + r)) 



-2(<5 2 (xK(x)vf(x + r)) 



(All) 



Performing as many angular integrals as possible, we are left with 
the following expressions, which we evaluate numerically. 

2 (<5, (x)v[(x + r)v r 2 (x + r)) = ^ J P(w)P(wy)w 3 j (wr)dwdy x 

-18y + 66y 3 + 66y 5 - 18y 7 + 9(-l + y 2 ) 4 7/ 
672y 3 

- {^"H 2 f J P(w)P(wy)w 3 J(^ n wr)dwdy X 
2y(9 - 109y 2 + 63y 4 - 27y 6 ) + 9(-l + y 2 ) 3 (l + 3y 2 )r 



672y 3 < A12) 
{6 2 (x)(v\(x + r)) 2 ) = -JL; J k 2 k 3 dk 2 dk 3 P(k 2 )P(k 3 ) x 

(-^Ll u (k 2 r)I u (k 3 r) + 

\ (| + ^j(4^ r I l2 (k 2 r)h, 2 (k 3 r) + (I - i/ £r )l 3 , (k 2 r)l 3 , a (k 3 r)) 

+ *(2tf r I u (k 2 r)I l , 3 (k 3 r) + I 3A (k 2 r)I 3A (k 3 r)XAl3) 



-2(5i(x)vf(x)v^(x + r)\ = ■ f k 3 dkdydxJ(jj. 2 tr ,kr)x 
x ' 2Stt J 

P(k ^+y 2 -2yx)P(ky) yX + ^ ~ 7X * 



1 + y 2 - 2yx 

.1.1.1 |..\ T- i ij - y.M.vu x yA)v\ 

2 



<A14) 



-2 (5i(x)v|(x)vf(x + r)) - 2 (<5 2 (x)vf (x)vf(x + r) 



— K f dyk 3 dkJ(n 2 ,kr)P(k)P(ky)x 
2ir J 

-2y(19 - 24y 2 + 9y 4 ) + 9(-l +y 2 ) 3 T 



<A15) 
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